Condensed Matter Physics, 200?, Vol. ?, No ?, pp. 0-?? 



U 









o 
o. 

^ Pressure-temperature phase diagram of 



^ : correlated hopping at half-filling 



o^. generalized Hubbard model with 



L.Didukh and V.Hankevych 

Ternopil State Technical University, Department of Physics, 

56 Rus'ka Str., UA-46001, Ternopil, Ukraine; E-mail: vaha@tu.edu.te.ua 



In the present paper pressure-temperature phase diagram of a generalized Hub- 
bard model with correlated hopping in a paramagnetic state at half-filling is 
determined by means of the generalized mean-field approximation in the Green 
Q . function technique. The constructed phase diagram describes metal-to-insulator 

Q ' transition with increasing temperature, and insulator-to-metal transition under 

the action of external pressure. The phase diagram can explain paramagnetic 
region of the phase diagrams of some transition metal compounds. 
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O . 1. Introduction 
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In the recent few years a generalized Hubbard model with correlated hopping 
has been used widely to describe strongly correlated electron systems (see papers [ 
|1|, 3 and references therein); the electron- hole asymmetry is a property of such a 
O ■ generalized Hubbard model as a result of the dependence of hopping integral on 

' occupation of the sites involved in the hopping process. Recently, this model has 

been extended to the case of a doubly orbitally degenerate band [ H . 
/\ • The generalized Hubbard model has much richer properties than the well-known 

j^ ■ Hubbard model [ ^, and an usage of the electron-hole asymmetry conception al- 

lows to interpret the peculiarities of physical properties of narrow-band materials 
which are not explained by the Hubbard model. In particular, the experimentally 
observed electron-hole asymmetries of metal oxides conductivity, of cohesive en- 
ergy of transition Sd-metaAs and of superconducting properties of high-temperature 
superconductors have been explained within the generalized Hubbard model with 
correlated hopping in papers [ §, §, 0) [01) [ S B respectively. 



Despite the fact that phase diagram of generalized Hubbard model has been stud- 
ied in works [ 0, [T^ |TB|, [T^, ^, jTB], [T^, researchers pay no attention to a determina- 
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tion of the model phase diagram in a paramagnetic state under the action of external 
influences, in particular pressure-temperature phase diagram. This task is related 
directly to the problem of metal-insulator transition description under the action 
of external pressure and temperature, namely the constructed pressure-temperature 
phase diagram of the model would allow to describe the observed metal-insulator 
transitions in narrow-band materials with changing of pressure and temperature. An 
interest to such transitions is caused by the theoretical point of view as well as the 
rich possibilities of its application (see, for example, monograph [ 1^ and review [ 



T^). Consequently, the goal of the present paper, being a continuation of previous 



work [ |2^ where the temperature-induced metal-insulator transition was studied, is 



to determine the pressure-temperature phase diagram of generalized Hubbard model 
with correlated hopping in a paramagnetic state at half-filling. On basis of this phase 
diagram we describe metal-insulator transition under the action of external pressure 
and temperature. 

2. Pressure-temperature phase diagram of the model 

Taking into account an external hydrostatic pressure p we write the model Hamil- 
tonian in the following form [ Q (in this connection see also Ref. [[2T[|): 



H = -ji^ af^ai^ + (1 + au)tY^ afaaja + ^2^ (af^ajarii^ + h.c. 

icF ija ija 

+Uj2ni^ni^ + -NVoKu^, (2.1) 

i ^ 

where i, j are the nearest-neighbours sites, [i is the chemical potential, a^, {aia) is 
the creation (destruction) operator of an electron of spin o [a =1; i) oii i-site (a 
denotes spin projection which is opposite to a), riifj = af^aia is the number operator 
of electrons of spin a on i-site, U is the intra-atomic Coulomb repulsion, t = t^ + Ti, 
with to being the matrix element of electron-ion interaction, Ti, T2 are the correlated 
hopping integrals (matrix elements of electron-electron interaction), the primes at 
the sums in Hamiltonian (2.1) signify that i ^ j- 

The latest term of the Hamiltonian has meaning of the elastic energy of a uni- 
formly deformed crystal, where k is the "initial" (purely lattice) bulk elasticity, A^ 
is the number of lattice sites, u = AVq/Vo is the relative change of the volume in 
unifrom strain (Vq is the initial unit-cell volume). Formulating the Hamiltonian we 



have used the result of paper [ ^: the dependence of a bandwidth W on relative 



change of the volume u in unifrom strain can be written in the form W = 2w{l + au), 
where w = z\t\ {z is the number of nearest neighbours to a site), a = ^^ < 0. 
We assume also that under the action of external pressure bandwidth changes only, 
and the matrix elements of electron-electron interaction (the correlated hopping in- 
tegrals and intra-atomic Coulomb repulsion) do not depend on relative change of 
the volume. 

As in papers [0, 20], using generalized mean-field approximation [ |^, 22 



an 



analog of the projection operation) in the Green function method we obtain for a 
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paramagnetic state at half-filling the energy gap width as 

AE = -(1 - 2d){w + w)[l + au] + -{Qi + Q2), (2.2) 

Qi = ^J[B{w -w){l + au) - Uf + \Uzt'{\ + (xu)^, (2.3) 

Q2 = ^\B[w - w))(l + au) + Vf + {Mzt'(\ + au)f, (2.4) 

where i? = 1 — 2(i + 4(i^, d is the concentration of polar states (holes or doublons) 
which has been calculated in Refs. [ |^, |2^, w = z\i\, i = t + 2T2, t' = t + T2; 
t and t are the terms describing hopping of quasiparticles within the lower and 
upper Hubbard bands (hopping of holes and doublons) respectively, t' describes 
a quasiparticle hopping between hole and doublon bands (the processes of paired 
creation and destruction of holes and doublons). 



According to the method proposed for the s(d)-f model in paper [ ^, the equi- 
librium value of relative change of the volume u is determined from the condition of 
minimum of the thermodynamic Gibbs' potential 

G = F + NpVo{l + u), (2.5) 

where F is the free energy. Using the known identity dF/du = (dH/du), Eq. ( |2.5| ) 
for the parameter u can be represented as 

_UiVpV„ = 0. (2.6) 



with H being Hamiltonian ( p.l|) . In the mean- field approximation passing to the 
space of quasi-momenta we get the following equation for relative change of the 
volume u: 

"w = -^^ Y. ^k(ai^^ak^) + rpVo, (2.7) 

where ai is the parameter which determines the quantity dW/dV, 2aiVo/W ~ 
0.1, r^0.05eV-i [0. 

Taking into consideration the fact that within the generalized mean-field approx- 
imation the first term of right-hand side of Eq. ( p.7| ) is equal to zero [ |2^ at the 



point of metal-insulator transition, we obtain the relation between relative change 
of the volume u and an external hydrostatic pressure p as 

au = rpVo. (2.8) 

Note that within the generalized Hartree-Fock approximation this equation is valid 
at the point of metal-insulator transition as well as in an insulating phase. 

To determine pressure-temperature phase diagram of the model we use for- 
mula (|2.2|) for the energy gap width and the expression for concentration of polar 



states calculated in Ref [ 20|. Let us consider, for instance, the Mott-Hubbard com- 



pound NiS2. This has two electrons half filling an €„ band, the half- width of initial 
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(uncorrelated) band of this crystal is wo 

-30 



z\to\ ^ 1.05 eV [ m, 13' and the initial 
unit-cell volume is Vq ~ 14.79 x 10"'^'^ m*^ [ ^. It shows the transition from the 
state of a paramagnetic insulator to the paramagnetic metal state at a hydrostatic 
external pressure of 46 kbar and room temperatures. Thus the transition occurs for 
a decrease in volume of about 0.4% with no change in crystal structure [ ^, 0. It 
also becomes metallic on alloying with Ni2Se, and the behaviour of this system is 
discussed later in this section. 

To calculate the model parameter U we fix one of the points {p = 22 kbar, 
T = 100 K) of the experimental curve in the phase diagram (the dashed- 
line curve of figure [l|) and find the value of intra-atomic Coulomb repulsion U 
at which the theoretical calculations within the present model reproduce this 
point. Thus, we obtain: ^ = -^ = 2.0168 for the correlated hopping parame- 
ters Ti = Ti/|to| = 0, r2 = T2/|to| = (these values of ti, T2 correspond 
to the Hubbard model), ^ = ,„^,f_^^, = 1.79107 for n = T2 = 0.1, and 



1.81437 at ri 



0, 



u 

u)o(l- 
^2 = 



0.1. Using these values of the model pa- 



LL = IL 

W Wq 

rameters we find the values of external hydrostatic pressure and temperature at 
which energy gap width is equal to zero (i.e. metal-insulator transition occurs). 
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Figure 1. 

The pressure-temperature phase diagram 
of metal-insulator transition determined 
within generalized Hubbard model with cor- 
related hopping for NiS2 in comparison with 
the experiment (the dashed curve): ti = 
T2 = (the upper curve); Ti = 0, r2 = 0.1 
(the middle curve); ri = r2 = 0.1 (the lower 
curve). PI (PM) denotes paramagnetic in- 
sulating (metallic) phase. 



The determined in this way 
pressure-temperature phase dia- 
gram of the model (figure [1|) de- 
scribes metal-insulator transitions 
in a paramagnetic state under 
the action of external pressure 
and temperature in NiS2, namely 
the constructed phase diagram de- 
scribes metal-to-insulator transi- 
tion with increasing temperature, 
and insulator-to-metal transition 
under the action of external pres- 
sure. Comparison of this theoret- 
ically determined phase diagram 
with the phase diagram of the com- 
pound NiS2 shows a good agree- 
ment between the theory and ex- 
periment. Besides, the theoretical 
calculations within the model re- 
produce the experimental data of 
paper [ ^ which point out a pres- 
ence of the energy gap width AE > 
in the ground state of NiS2 and 



in absence of an external pressure. 
The phase diagram shows that taking into account correlated hopping allows much 
better description of these experimental data than the Hubbard model; this testifies 
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also much better physics of the present model and the important role of correlated 
hopping. 

Analogous phase diagrams can be constructed for the other compounds: 
(Vi_,Cr,)203 [ 0, H, NiS2-.Se, [ H, and Yi.^Ca.TiOa [ IT], H exhibiting 
such metal-insulator transitions. For example, the materials (Vo.96Cro.o4)203 shows 
a metal-insulator transition at a hydrostatic external pressure of 13 kbar and room 
temperatures; the transition occurs for a small decrease in volume of about 1% with 



no change in crystal structure [ 25]. In (Vi_a;Mi,.)203 (with M = Cr, Ti) the addition 
of Ti^"*" ions to V2O3 leads to insulator-to- metal transition, whereas the addition of 
Qj^3+ JQj-^g j-gg^i^g in inetal-to-insulator transition. The simplest explanation is that 
the substitution of V^^ ion for Cr^"*" ion leads to a band narrowing; the Cr^"*" ion 
is a localized impurity and it deletes a state from the 3(i-bands, deleting a state is 



equivalent to a band narrowing or an external pressure decreasing [ |18[ which drives 
the system towards the insulating phase. Likewise the addition of Ti^"*" impurities is 
equivalent to an external pressure increasing. 

In the Mott-Hubbard compound NiS2-xSe, electron hoppings between the sites 
of Ni occur by the chalcogenide sites (this is caused by peculiarities of the pyrite 
crystal structure [ 0), the substitution of S^~ ion for Se^~ ion in NiS2 leads to an 
increase of wave functions overlapping, consequently the probability of an electron 
hopping increases which is equivalent to a band broadening or an external pressure 
increasing. Therefore, the pressure-temperature phase diagram constructed for NiS2 
can describe the experimental composition-temperature phase diagram [ ^, ^ of 
the compound NiS2-xSe2 



-'.• 



Note also that to construct phase diagram of the system Yo.eiCao.agTiOs we 
have to generalize the previous results obtained at half-filling to the case of a non 



half-filled band because this compound is characterized by such a band [ ^ . 

In conclusion, in the present paper pressure-temperature phase diagram of the 
generalized Hubbard model with correlated hopping in a paramagnetic state at 
half-filling has been determined. The constructed phase diagram describes metal-to- 
insulator transition with increasing temperature, and insulator-to-metal transition 
under the action of external pressure. Comparison of this theoretically determined 
phase diagram with experimental data, in particular with the phase diagram of the 
compound NiS2 shows a good agreement between the theory and experiment. We 
have found that taking into account correlated hopping allows much better descrip- 
tion of these experimental data than the Hubbard model; this testifies also much 
better physics of the present model and the important role of correlated hopping. 

The determined pressure-temperarure phase diagram of the model can explain 
paramagnetic region of the phase diagrams of the transition metal compounds: the 
systems NiS2-a;Sea; and (Vi_a;Cra,)203, calcium doped YTiOa. 
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